Negative differential resistance in scanning tunneling microscopy: simulations on 

Cgo-based molecular overlayers 



Frederico D. Novaes 1 , Manuel Cobian 1 ' 2 , Alberto Garcia 1 , Pablo Ordcjon 3 , Hiromu Ucba 4 , and Nicolas Lorente 3 
1 Instituto de Ciencia de Materiales de Barcelona (CSIC), E-08193 Bellaterra, Spain 
2 LPMCN, UCBL, CNRS, UMR5586, 69622 Villeurbanne, France 
3 Centre d'Investigacio en Nanociencia i Nanotecnologia (CSIC-ICN), E-08193 Bellaterra, Spain 
4 Division of Nano Science, Graduate School of Science and Engineering, University of Toyama, Toyama, Japan 

(Dated: January 20, 2011) 

We determine the conditions in which negative differential resistance (NDR) appears in the C60- 
based molecular device of [Phys. Rev. Lett. 100, 036807 (2008)] by means of ab-initio electron- 
transport simulations. Our calculations grant access to bias-dependent intrinsic properties of the 
molecular device, such as electronic levels and their partial widths. We show that these quantities 
depend on the molecule-molecule and molecule-electrode interactions of the device. Hence, NDR 
can be tuned by modifying the bias behavior of levels and widths using both types of interactions. 



Since the creation of the first tunnel diode a non- 
linear I—V characteristic, particularly showing negative 
differential resistance (NDR) is the cornerstone in two- 
terminal devices [H, Q ■ NDR permits the electronic cur- 
rent of these devices to be modified by the applied bias 
in a complex way. This active control of the current mo- 
tivates a long search for NDR in molecular devices [3-flo| 
where the scanning tunneling microscope (STM) plays 
an important role [H, @] . The STM can vary applied bias 
and tunneling current in an independent manner, yield 
atomically precise data and still be a two-terminal de- 
vice. More recently, STM studies of Cgo molecules have 
spurred a lot of interest because of the NDR capabilities 
of Cgo-based devices 

As early as the first NDR cases were obtained in STM 
junctions, the existence of sharp resonances on both elec- 
trodes were used to explain the decrease of current as 
the bias increased @, 0, 14 1. However, Grobis and co- 



workers |ll| explained their own NDR data by uncov- 
ering yet another mechanism: the voltage-dependent in- 
crease of the tunneling barrier height. Recently, one more 
mechanism has been advanced in order to explain the de- 
pendence of NDR on the tip's material, namely orbital 
matching between molecule and tip (lo| . This last mech- 
anism has been thoroughly studied and validated, while 
showing the need of sharp STM tips [ltij . However, in 
most molecular devices, the accurate evaluation of the 
voltage drop across the device is crucial. As a matter 
of fact, Tu and co-workers [13] showed the paramount 
importance of the actual voltage-dependent barrier for 
the presence or absence of NDR. This spread of mech- 
anisms shows that many parameters change among dif- 
ferent molecular junctions and it is of fundamental im- 
portance to understand the key ingredients in NDR and 
thus the final functionality of possible molecular devices. 

In this letter, we study the molecular device of Ref. (l2| 
unrave ling the mechanisms leading to the measured 



(TPA) molecules. Despite the complexity of the sim- 
ulations, we can extract the molecular levels and their 
coupling with the electrodes while describing their bias 
dependence. This allows us to have unprecedented in- 
sight in the actual physics of a molecular device with a 
quantitative description. We show that we can tune the 
size of the NDR signal by modulating the interplay of 
the bias dependence of the levels with the bias depen- 
dence of their partial widths which ultimately depend on 
molecule-molecule and molecule-electrode interactions. 



In the Landauer formalism (18(, the electronic current, 
I, is given in terms of the electron transmission function 
T(E, V) with E the electron energy and V the applied 
bias (e and h being the electron charge and Planck's con- 
stant, respectively), by: 
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J dE [f R (E, V) - f L (E, V)} T(E, V), 



NDR 12]. The simulated molecular device consists 
of Ceo molecules partially decoupled from a Au (111) 
substrate by spectator 1,3,5,7-tetraphenyladamantane 



(1) 

where the Fermi filling factors of the left and right elec- 
trodes, Jl and fn include the bias, and correspond to 
the chemical potential of the left and right electrodes 
asymptotically inside each electrode. NDR corresponds 
to negative conductance. From Eq. ([I]), we can derivate 
with respect to the bias to obtain the conductance. We 
obtain two terms, one comes from the derivative of the 
Fermi filling factors and the other one corresponds to 
the derivative of the transmission. The derivative of the 
Fermi factors is a positive number, hence, only the deriva- 
tive of the transmission gives rise to NDR. It is then, the 
behavior of the transmission function with bias that de- 
termines the appearance of NDR. 

A quantitative way to compute Eq. ([I} is to use non- 
equilibrium Green's functions (NEGF) self-consistently 
solved with Poisson's equations and fixed chemical po- 
tentials in the asymptotic region inside the electrodes. 
This has been done using the Transiesta package (l9| . 
where the electronic structure for the NEGF is evalu- 
ated using density functional theory (DFT) in the local 
density approximation (LDA). To optimize our calcula- 
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FIG. 1. Geometry and voltage drop for a 0.3-V bias in the 
four systems analyzed in this work: (a) a Ceo-TPA overlayer 
on a semi-infinite Au(lll) an a second semi-infinite Au(lll) 
representing a flat STM tip, the tip- molecule distance is 6.7 
A, (b) the same overlayer under a pyramidal tip formed by 
4 atoms and a semi-infinite Au(lll) surface, the tip-molecule 
distance is 6.6 A, the tip apex is off the center of the topmost 
Ceo hexagon, (c) same as (b) but a capped pyramide with 
three atoms, the tip-molecule distance is 6.9 A, (c) same as 
(a) plus a gold adatom, the tip-molecule distance is 6.9 A, the 
tip apex is off the center of the topmost Ceo hexagon. The 
voltage drop is shown on two planes cutting through the Ceo 
and the TPA molecule respectively. The plotted isolines are 
0.0083 V apart. 



tions, we use a special surface basis-set that describe the 
electronic properties at the metal- vacuum interface [2l| . 
Hence, the present electronic structure of the molecu- 
lar overlayer on Au (111) reproduces the LDA results 
of Ref. [12J. The size of the electronic problem is very 
large since 5 active layers are used for right electrode, R, 
and four for the left one, L. Each layer contains 28 Au 
atoms [T3|. The junctions include one Cgo molecule and 
one TPA, plus the STM tip. Figure Q] shows the geometry 
and the voltage drop at the STM junctions considered in 
this work: (a) a flat tip, (b) a pyramide, (c) a capped 
one, and (d) an adatom. The atom-terminated tips are 
placed in a low-symmetry point over the Ceo molecule, 
in order to avoid an unrealistic high-symmetry configu- 
ration. The two topmost layers of the R electrode and 
the topmost one of the L electrode are relaxed, the other 
three are fixed to the bulk position, and a recursion al- 
gorithm is used to reproduce semi-infinite electrodes (l9j 
while allowing for the correct electric field screening and 
consequent Friedel oscillations. 

Figure [T] also shows the voltage drop across the tun- 
neling junction when the external bias is 0.3 V. We see 



that most of the bias drop takes place in the vacuum re- 
gion. At this bias, the LUMO [20j of Ceo is close to the 
asymptotic chemical potential of the tip. Hence, Ceo is 
largely conducting and the bias drop is smaller than in 
the vacuum gap: the voltage drops a tenth of the total 
bias just between the top of the molecule and the metal- 
lic substrate. The molecular states are then subjected 
to shift under the action of the bias. However, its part- 
ner TPA molecule has a large HOMO-LUMO gap which 
makes it basically inert. The voltage drop takes place 
across the molecule and it amounts to a seventh of the 
total bias. Hence, TPA behaves as the dielectric mate- 
rial of the molecular device that changes the charging 
energy of the active molecule, Cgo- Calculations with 
and without TPA show that the charging energy of Ceo 
is reduced by 0.2 eV in the presence of TPA. Finally, as 
the tip becomes sharper, the electric field gets more lo- 
calized. As a consequence, the electric field is large close 
to the tip and the voltage drop is faster in the vacuum 
gap while smaller at the molecular site. Hence, by sharp- 
ening the tip we obtain an effect similar to increasing the 
vacuum gap which enhances the tunneling-barricr bias 
dependence (l7j . 

The transmission behavior with applied bias varies 
rapidly as the tip changes. Figure [2] shows the electron 
transmission T of Eq. ([1]) as a function of the electron en- 
ergy E, as the bias between tip and surface is increased. 
For a fixed bias, V, the transmission function can be re- 
produced by a Breit-Wigner (22j formula: 
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The full-line curves in Fig. [2] are fits to Eq. ([2]) and dots 
are the Transiesta-evaluatcd transmissions. The agree- 
ment is excellent permitting us to extract the partial 
widths, Tr and Tl as well as the LUMO level Eq as 
a function of applied bias, see Fig. [3J 

The bias dependence of Eq, Fig. [3J shows that the 
LUMO level is not pinned to the substrate's electronic 
structure because the bias is not entirely droping in the 
vacuum gap. This is due to the molecule-molecule in- 
teractions of Ceo with surrounding TPA, that decouple 
from the substrate and polarize Ceo- The fiat-tip causes 
the largest Eq shift both at zero-bias and at finite bias, 
showing that despite the 6.7-A vacuum gap, there is in- 
teraction between molecule and tip. The tip also modifies 
the coupling of the molecule with the substrate. Indeed, 
the partial width due to the substrate, T^, can be divided 
in two values, Fig. [3J one corresponding to the blunt tips 
(a) and (c) and the other one to the atomic-like tips (b) 
and (c). Due to the stronger coupling of the molecule 
with the Au(lll) substrate, Tr is larger and has a much 
smaller U-dependence than Tl that is the coupling with 
the tip. The behavior of with V can be basically 
understood by the shift of the LUMO level, Eq, As the 
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FIG. 2. Electron transmission as a function of the electron 
energy (E) the same junctions of Fig. [1] (a) the fiat-tip junc- 
tion (b) the pyramidal one (c) the capped one and (d) the 
adatom one. Filled circles are Transiesta results, and the 
full line is the fit to Eq. (J2j) . The peaks shift with increasing 
bias, hence the lower-energy peak corresponds to 0.0 V, the 
next one to 0.4 V, 0.8 V, 1.2 V and 1.6 V respectively. Please, 
notice the different transmission scales on both sides. 



bias increases, Eq shifts closer to the vacuum level of the 
substrate, slightly increasing the partial width r^. How- 
ever, is strongly dependent on the actual density of 
states (DOS) of the tip's electronic structure. Sharp tips, 
(b) and (d) , show a non-monotonic behavior of Tl with 
V. This is also seen in the electron transmission since 
by virtue of Eq. ©, the transmission at Eq is roughly 
proportional to T^. 

Figure [2] (a) is the transmission function for the flat 
tip case. We see that the transmission maximum, cor- 
responding to the LUMO, shifts to higher energies as 
the bias is increased. This is due to the partial bias 
drop at the center of the molecule and is reproducing 
the behavior described in Ref. [ItJ- More interestingly, 
the maximum height drops with bias. This is the behav- 
ior described in Ref. [TTJ] causing NDR. And indeed, as 
can be seen in Fig. 2] (a) NDR is found for this system. 
This fast decay of the transmission is due to the increase 
of the tunneling barrier with increasing bias 
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For the 4-atom pyramidc tip, Fig. [5] (6), the maximum 
of the transmission peaks when the bias is at ^0.8 V. 
This is totally different from the behavior of Fig. [2] (a) 
and the one described in Ref. [TTJ] . The maximum of the 
transmission with bias can be traced back to a resonance 
of the 4-atom pyramidal tip of atomic origin. This is 
found at -0.45 eV in the DOS of the tip. Then, the 
strong variation of Tl with bias, Fig [3l is due to the 
resonance in the DOS. When the tip is capped, Fig. Q] 
(c), the corresponding recovers a monotonic decrease 
with bias as for the flat-tip case (a), Fig. [31 Hence, the 
transmission is very similar to the case of the flat-tip. 
Finally, case (d) where one adatom is added to the flat 
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FIG. 3. Bias dependence of the LUMO level Eo with respect 
to the substrate's Fermi energy, its partial width due to its 
coupling to the substrate, Tr and the partial width due to 
the coupling to the tip, Tl- Symbols: diamonds correspond 
to the flat tip of Fig. [1] (a), triangles to the 4-atom tip, (b), 
squares to the capped tip, (c), and dots to the adatom tip, 
(d). 



tip shows the same kind of non-monotonic behavior in 
as in (b) and a clear maximum in both, and the 
transmission peaks for the LUMO around 0.8 V. 

Figure 2] presents the I-V characteristics for the four 
systems of Fig. Q] The NDR features are large for the 
four-atom pyramidal tip, Fig. [4] (b) and small for the 
flat tip, Fig. |4] (a). Hence, the largest NDR is obtained 
when the partial widths present a resonant behavior with 
the applied bias, as the one due to electronic resonances 
in the electrodes, tip (b). The monotonic bias depen- 
dence of partial widths due to the barrier dependence 
on the bias, case (a), is slower leading to weak NDR. It 
is interesting to study the capped pyramidal tip, Fig. @] 
(c), where NDR is absent and the I-V curve is closer to 
the flat-tip one. Here, the partial width dependence is 
slow, Fig. [3J with increasing bias because the tip's DOS 
is smoothly increasing. Despite the fact that the barrier 
height increases with bias, reducing the partial width, the 
increase of DOS is large enough to partially compensate 
the larger barrier and the transmission only drops slowly, 
Fig. [2] (c). NDR disappears. However, if one adatom is 
added to the flat tip, a bigger NDR effect is obtained be- 
cause the sharper tip DOS forces a resonant-like Tl and 
hence a peak-shaped transmission. It is then this last 
effect, rather than the enhancing of the bias dependence 
of the tunneling barrier, Fig.Q] that drives NDR in sharp 
tips. 

The shape of the I-V characteristics changes with tip. 
While the flat tip has a sharp onset and a sharp peak 
at 0.6 V, the atom-terminated tips have a broad peak 
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FIG. 4. I-V characteristics for the four cases studied here in 
the same order as Fig. [T] Dots are the Transiesta results, 
and full line the results of Eq. fl]) using the Breit-Wigner fit, 
Eq. J2j), interpolating the LUMO and LUMO+1 levels and 
partial widths for bias in between the Transiesta calcula- 
tions. The fit excellently reproduces the LUMO onset in the 
current, however errors accumulate in between the LUMO 
and LUMO+1 onsets where the Breit-Wigner fit is worse. 
The current peak-to-valley ratios (dots) are (a) 1.07, (6) 1.4, 
(c) - , (d) 1.1, the experimental one [lJl is 1-3. 



centered at 0.8 - 0.9 V with a ~ 0.5- V width. The ex- 
perimental one is a broad peak centered at 1.1 V with 
a ~ 0.5-V width [l2j]. This fact and the current peak- 
to-valley ratio, Fig. 21 render the pyramidal tip in very 
good agreement with the experiment. Moreover, the ex- 
perimental distance between the LUMO and LUMO+1 
onsets [H is 1.5 V and the one in Fig. H is 1.2 V. The 
actual I-V shape is then determined by the bias depen- 
dence of the level, Eo, and of the partial widths, here 
Fl. Indeed, the shape of, for instance, the LUMO onset 



23j, while the NDR 



is given by the U-dcpendcncc of Eq 
drop is dominated by the U-depcndence of Tl ■ 

In summary, we have performed NEGF calculations for 
the experimental setup of Ref. [l2| • Our simulations have 
permitted us to rationalize the appearance of NDR in this 
system in terms of molecular levels and partial widths be- 
cause we have been able to extract their bias dependence. 
We have shown that NDR can be tuned by modifying 
molecule-molecule and molecule-electrode interactions. 
The molecule-molecule interactions induce a reduction 
of the overwhelming presence of the substrate permit- 



dependence. On the other hand, molecule-electrode in- 
teractions determine the molecular partial widths. NDR 
effects are maximum when the partial widths vary rapidly 
with the applied bias. This dependence is largest in the 
presence of electronic resonances on the electrodes such 
as the ones obtained for sharp STM tips. Finally, our re- 
sults show that the overall shape of the I-V characteristic 
contains relevant information on the molecule-molecule 
and molecule-electrode interactions. 
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